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ABSTRACT 



The Y.T. Lee Array for Microwave Background Anisotropy (AMiBA) started scientific operation in 
early 2007. This work describes the optimization of the system performance for the measurements of 
the Sunyaev-Zel'dovich effect for six massive galaxy clusters at redshifts 0.09 — 0.32. We achieved a 
point source sensitivity of 63 ± 7 mJy with the seven 0.6 m dishes in 1 hour of on-source integration in 
2-patch differencing observations. We measured and compensated for the delays between the antennas 
of our platform- mounted interferometer. Beam switching was used to cancel instrumental instabilities 
and ground pick up. Total power and phase stability were good on time scales of hours, and the 
system was shown to integrate down on equivalent timescales of 300 hours per baseline/correlation, 
or about 10 hours for the entire array. While the broadband correlator leads to good sensitivity, the 
small number of lags in the correlator resulted in poorly measured bandpass response. We corrected 
for this by using external calibrators (Jupiter and Saturn). Using Jupiter as the flux standard, we 
measured the disk brightness temperature of Saturn to be 149;';^2 K. 

Subject headings: cosmic microwave background — galaxies: clusters: general — instrumentation: 
interferometers 



1. INTRODUCTION 

The angular power spectrum of cosmic microwave 
background (CMB) anisotropics carries a wealth of in- 
formation on the physical processes in early epochs of 
the universe. A comparison of theoretical models with 
accurate measurements of CMB anisotropics thus con- 
strains the fundamental cosmological parameters and 
models for cosmic structure formation. On larger angu- 
lar scales, the temperature anisotropics are dominated by 
primary CMB fluctuations, whereas on smaller angular 
scales secondary effects such as the Sunyaev-Zel'dovich 
(SZ) effects due to galaxy clusters dominate over pri- 
mordial anisotropics. The amplitude and location of 
the peak in the thermal SZ power spectrum are par- 
ticularly sensitive to the amplitude of the primordial 
matter power spectrum, represented by the normaliza- 
tion (Tg: as well as the thermal history of the hot intra- 
cluster medium. Th e Cosmic Background Imager (CBI, 
iPearson et al.ll2003[ ) and Arc minute Cosmolog y Bolome- 
ter Array Receiver (ACBAR, iKuo et al.ll2004l ) measured 
the CMB temperature power spectrum at large angular 
multipoles of / ^ 3000. While the CBI detected an excess 
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power over the theoretical prediction from the standard 
cosmological model, the ACBAR result has a larger error 
bar and is consistent with both an excess and no excess. 
To date the uncertainties of the high-Z measurements re- 
main large. More accurate measurements on large angu- 
lar scales around and beyond I = 30 00 are required to 
better constrain the value of gs (e.g . , iBond et al] l2005t 
iGoldstein et"aI1[2f)0a iLin et al.l[20f)l . 

The Y.T. Lee Ar ray for Microwave Background 
Anisotropy (AMiB A, iHo et al.l l2009l : IChen et all l2009l : 
iKoch et all l2009af) is designed to measure CMB 
anisotropics on these multipole scales. The AMiBA is 
located on the volcanic mountain Mauna Loa, Hawaii, at 
an altitude of 3400m. The array observes with a single 
sideband in 86 — 102 GHz, or at roughly 3 mm wave- 
length, with cooled HEMT low noise amplifiers (LNA). 
Each of the seven receivers measures two linear polar- 
izations (X and Y) and produces two corresponding IF 
channels (each 2—18 GHz). Out of the four possible 
cross-correlations with a pair of receivers, AMiBA em- 
ploys a switching system to form either the (XX*, YY*) 
or the (XY*, YX*) product at the same time. Note that a 
circular polarizer is being developed so that AMiBA can 
choose to measure either the (LL*, RR*) or the (LR*, 
RL*) cross-correlations in the future. There are thus 
21 baselines and 42 instantaneous correlations for the 
seven-element array. The correlation is further divided 
into complex visibilities in two frequency bands using an 
analog four-lag correlator ()Li et al.ll2004l ). 

All antennas and receivers are mounted on a 6 m plat- 
form so that antennas can be closely packed without 
issues with shadowing and collision. In the 2007 and 
2008 seasons, observations were made with 60 cm diam- 
eter dishes close-packed in the center of the platform. 
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IWu et all (|2009D present details of the observations and 
analysis of six massive clusters. 

In this paper we describe how the system performance 
was optimized for these targeted observati ons. Two com- 
panion papers discuss the data integrity (Nis hioka et al.l 
[2009() and the C MB and foregrou nd uncertainty in the SZ 
flux estimation (|Liu et al.ll2009f ). Combined with pub- 
lished X-ray parameters, the SZ fluxes of six clusters 
were u sed to measure the Hubble paramete r (iKoch et al. 
200%") and to examine the scaling relations (|Huang et al.l 
2009,). Subaru weak lensing data for four of the clusters 
were analyzed with the SZ measure ments to derive the 
baryon fraction (jUmetsu et al.l[2009^ ). 

This paper is organized as follows. Critical issues such 
as the noise temperatures, delay corrections, stability, 
spurious signal removal and characteristics of the cor- 
relators are described in fJH discusses the losses of 
the system, the calibration errors, and the integration of 
noise. Finally 21 summarizes our conclusions. 

2. OPTIMIZING INTERFEROMETER PERFORMANCE 

Prior to and during the 2007 observing season, commis- 
sioning activities identifi ed parts of the o perations which 
needed to be impro ved (|Lin et al.l 120091 ). In particular, 
iHuang et al.l ([2 00 8*) reports on the deformation of the 
platform which can affect the performance of the inter- 
ferometer. Fortunately, these platform errors are repeat- 
able and can be modeled. Their effects on poi nting, radio 
alignm ent, and phase errors are discussed in iKoch et al.l 
i2009a). For AMiBA operations in 2007-8 these effects 
were minimal. In this paper, we concentrate on other 
areas of the interferometer performance which were op- 
timized. 

2.1. System Temperature 

To understand the gain stability of AMiBA, we first 
measured the receiver stabilities. The system temper- 
ature is monitored by a set of sky-dips in total power 
mode. The total power output from each IF channel can 
be approximated by 

PiF — gkB[Trx + Tdish + Tcmb 

+ Tatm/ sin{el) + Tgndiaz, el)], (1) 

where g is the power gain, k is the Boltzmann constant, B 
is the bandwidth of each IF channel, and the T's denote 
the noise temperatures from the receiver (rx), antenna 
(dish), CMB (cmb), the atmosphere (atm), and ground 
pickup {gnd}. A hot /cold load measurement is used to 
calibrate gB and Trx- The receiver noise temperatures 
are 55- 75 K ( Chen ct al. 2009). Fitting the total power 
to P = Pq + Pi/ sin(e^) lumps the contributions into sky- 
like (Pi) and receiver-like (Pq) parts plus some residual 
contributions from the ground. The measurements show 
that the total receiver-like noise temperature is about 
1(7 ~ 5 K higher than Trx- The sky- like part is approxi- 
mately 15 K at zenith in typical observing conditions. In- 
cluding Tcmb, the system temperatures away from zenith 
are about 80 — 100 K. Repeated hot/cold load measure- 
ments of the receiver noise temperatures show that T^x 
is stable within the measurement error (~ 5 K). Hence, 
by monitoring the system temperature using sky-dips, 
we can reject inferior sky conditions and unstable instru- 
ment behavior. 



2.2. Delay Correction 

Since AMiBA is a coplanar array there is no fringe 
rotation in a tracking observation. Fringes occur when 
a source moves across the field of view (fov) creating a 
geometric delay. The fov of AM iBA equipped with 0.6- 
m dishes is 23' (|Wu et al.l I2009D . The requirement on 
delay trimming is that the source delay should remain 
within the sampling range of the lag-correlator, which is 
±50 ps. As the source delay approaches the limit of sam- 
pling range, the error in the recovered visibility becomes 
larger with a consequent rapid drop in sensitivity. To 
allow a 2-m baseline to observe a 23' fov, which corre- 
sponds to a delay range of ^ ±22 ps, the instrumental 
delay was specified to a tolerance of ±20 ps. 

To measure the delay for each correlation, all dishes 
were removed and a noise source was mounted between 
receivers (e.g. Anti and Ant2). A fringe is generated 
when the noise source moves from Anti toward Ant2, 
simulating a fringe due to a celestial source. 

L{X, Ta) = ^ d/R'(/)e~^2.[(/+/,o)^+/(r.-ri+r„)]^ 

. . (2) 

where x is the displacement of noise source, / is the IF 
frequency, and R' is the complex response function of the 
baseline excluding the linear part of the phase due to lags 
(tq, a = 1...4) in the correlator. TZ takes the real part of 
the expression and is done implicitly whenever necessary 
hereafter, ri and T2 represent the instrumental delays in 
the IF's of Anti and Ant2. The fringe envelope peaks 
when ^ = Ti — T2 — Ta- The relative delay ri — T2 is 
measured with respect to the central lag (with Ta = 0). 
Equation ^ is usually referred to as the lag output or 
the lag data throughout this work. 

Wc found the instrumental delays for all IF channels 
using relative delay measurements. Short cables were 
then inserted into each IF for compensation. After this 
trimming procedure the residual delays were measured 
by fitting fringes for the Sun without the dishes, model- 
ing the fringes as the convolution of the observed point 
source fringe with a circular disk. The differences be- 
tween observation and model are consistent with resid- 
ual delays of ±15 ps (RMS). Except for the delays due to 
platform deformation, the delays between antennas were 
therefore well controlled. 

2.3. Bandpass Shape Measurement 

The AMiBA correlator has four lags and outputs two 
spectral bands to cover the 2 — 18 GHz band. Know- 
ing the bandpass shape is an important aspect of ob- 
taining good visibilities using this type of correlator (see 
next section Because the analog correlator con- 

tributes significantly to the bandpass shape, we adopted 
a baseline-based measurement approach. The Fourier 
transform of the fringe L{x,Ta) against x is used to de- 
termine R' for each baseline, with a spectral sampling 
of about 0.8 GHz. Each lag output is transformed inde- 
pendently. Fig. [1] displays the gain and phase responses 
of all valid measurements after the four lag outputs are 
averaged together. Averaging the phase responses of the 
four lags, the delays are canceled leaving only the com- 
mon mode of the spectral variation. The gain responses 
of the four lags are summed and then normalized such 
that the averaged gain between 2—18 GHz is set to 1. 
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The conversion from observed fringe rate to the RF 
frequency is proportional to the noise source translation 
speed. We believe a ±1 % jitter is present in the trans- 
lation stage we used, which introduced roughly ±1 GHz 
uncertainty in the response frequency. This causes one 
of the major problem in developing an accurate visibil- 
ity extraction method from our lag data (see ^2.4p . An 
improved measurement setup involving simultaneous in- 
jection of a single frequency source is being developed to 
achieve a higher accuracy. 

The effective band widths, defined as B = 
|/d/R'p//ci/|R'p, are insensitive to the uncer- 
tainty in the response frequency. Based on our bandpass 
measurements, the effective bandwidths of the AMiBA 
correlators are calculated and shown in FiglS] They 
generally fall in the range of 7 — 13 GHz. 

2.4. Extracting the Interferometer Visibilities 

Several approaches can be used to convert the four 
measured lags of the AMiBA correlator into complex vis- 
ibilities in two bands over the 16 GHz bandwidth. We 
find that the inaccuracies inherent in this inversion need 
to be corrected by externa l calib ratio n. Here we ad opt 
the formalism of IWu et all (|2009[ ) fsee lLi et al.ll200 4 for 
an alternative formalism). The lag output in equation ^ 
can be expressed in matrix form as La = KafcV^'"'^, where 
YfT'^ is the source visibility. Subscript a indexes the 
^lag = 4 lags, and subscript k indexes the Nf discretized 
frequency samples fk, where Nk is usually much larger 
than Niag. 

The transformation relies on a kernel Kak, which is an 
estimate of the response matrix Mafe- The kernel is inte- 
grated in frequency into two bands Kac, where c = 1...4 
indexes the real and imaginary parts of the two bands. 
We use the inverse of the integrated kernel to construct 

the raw visibihty V™"' = K'^a^a- 

Ideally we would like Kak = ^ak so that V™"" is clos- 
est to V^'''^. However, when Kak is an inaccurate repre- 
sentation of Rq^:, due to measurement errors, variations 
with temperature or time, insufficient spectral resolution 
in the measurement, or insufhcient information about 
the response, errors in visibilities occur. Fig. [3] demon- 
strates the calculation of raw visibility using simulated 
drift scans in three cases when (1) the kernel is the ex- 
act response, (2) there are measurement errors, and (3) 
there is no knowledge about the response. The correct 
visibility should appear as a Gaussian in amplitude with 
a linearly increasing phase. It can be seen that case (1) 
recovers the correct result, whereas deviations from this 
form increase with decreasing accuracy of the kernel. We 
therefore must obtain a calibrated visibility from the raw 
visibility Vg°' = CbcV™™, where b has the same index 
range as c, and Cbc is the calibration matrix, which can 
be obtained by comparing the raw visibility of a planet 
(the calibrator) to the theoretical visibility. 

In the analysis of data taken in 2007 and 2008, the 
flat kernel (right-most function in Fig. [3]) was assumed, 
and planet calibrations were applied. We have estimated 
the errors introduced by external calibration by running 
simulations on point source models. This error is on 
the order of ±2% (la) in the absolute fluxes with no 
detectable bias. This is small compared to the thermal 
noise and the measurement errors on the planet itself. 



2.5. Stability 

The stability of the system was examined by measuring 
the variation in visibilities for a few bright planets during 
local times 8 pm to 8 am, as normally used for observing. 
For this test, the ephemerides of the planets were taken 
at the beginning of each track but not updated during 
the observation. This causes a pointing error that in- 
creases to about an arcminute over 12 hours. To account 
for this, two sets of visibility data for each planet were 
chosen as calibrating events. A linear interpolation was 
used to remove the linear drift. For data without brack- 
eting events, the nearest calibration was used. Fig. [5] 
shows an example of a stability measurement. The gain 
stability was found to have an RMS variation around 5%, 
and the phase to have an RMS variation around 0.1 rad. 
The measurements also reveal that the phase response is 
more sensitive to changes in environment than the gain 
response, especially in the first hour after shelter open- 
ing. 

FiglS] plots the flux of Jupiter recovered from the data 
set used in FigH) Data was calibrated by the first mea- 
surement at UT 12h (not plotted). The recovered flux 
varied within ±4% of the expected flux until sunrise. 
Calculation of the calibrator flux is discussed in ^3.2[ 

Based on the stability measurements, we chose to use a 
calibration interval of two to three hours, to give calibra- 
tions good to about 5% in gain and 0.1 rad in phase for 
each baseline. Calibration requires ^ 10% of telescope 
observing time. 

2.6. Minimizing Instrumental and Ground Pickup 

The signal in the lag output should be constant when 
AMiBA tracks a source. However, the weak signal we 
measure is susceptible to slowly-varing contamination. 
The system is designed with a phase switching and de- 
modulation scheme to remove contamination such as a 
common mode leakage in the IF paths. To modulate the 
signal, we use a PIN switch to change the LO between 
two carefully adjusted delay lines, and thus changes the 
phase of IF signal by 180°. The demodulation is done 
in the readout process. The aim was to remove con- 
tamination between the down-conversion mixer and the 
correlator readout. 

However, mixers in the correlator can pick up higher- 
order signals such as |£'ip|£^2p in addition to their nomi- 
nal output, which is proportional to E1E2 or \E\'^ ^ where 
Ei stands for the voltage from Ant^. If the power of the 
IF signal is modulated by the phase switching pattern, 
then this higher order response can generate an output 
that is coherent with the demodulation pattern and be- 
comes a spurious signal. This effect is, indeed, seen in 
AMiBA, where phase switching of the LO can result in a 
power difference as large as 0.3 dB. The LO power mod- 
ulation is carried through to the IF in varying amounts 
depending on the mean LO power level at the mixer. The 
IF modulation can be undetectable for optimally-tuned 
mixers, but is up to 3 dB for under-pumped mixers. 

To reduce the IF modulation, the LO drive level is 
optimized for minimum conversion loss for each mixer. 
Modulation of the LO power would then have the least 
impact on the IF power. Furthermore, to reduce the 
drifting of LO level with ambient temperature changes, 
the final amplifier and frequency doubler in the LO chain 
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Fig. 1. — The complex responses of AMiBA. Responses include effects from the RF components, IF components, and the analog correlator. 
Top and bottom panels display the gain and phase responses respectively. Each line represents one cross-correlation (XX on the left and 
YY on the right) of a pair of receivers. Vertical dashed lines indicate RF frequencies of 84 and 102 GHz (IF frequencies and 18 GHz). 
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Fig. 2. — Effective bandwidths of the AMiBA correlators calcu- 
lated from the bandpasses displayed in Fig[l] The percentages are 
based on a nominal input bandwidth of 16 GHz. 

are operated in the soft saturation regime. Additional 
protection is provided by temperature controls installed 
before the 2009 observing season. 

Spurious signals external to the system, such as ground 
pickup, will still affect the data. We used a subtraction 
scheme similar to the one used by CBI (iPadin et al.ll200a 
to suppress the slowly-varying signals. In practice, we 
have found that the spurious signal in individual patches 
in ~5 hrs integration can be as high as ±7 Jy/beam, but 
that after subtraction a cluster with brightness 0.3 Jy can 
be detected at 9cr level in 11 hours (i.e., with 5.5 hours 
on-source integration). T he observing str ategy and the 
data analysis are given in iWu et all (pool . 

3. ACHIEVED SYSTEM PERFORMANCE 



3.1. Overall Efficiency 

For each baseline, losses include the antenna loss, an- 
tenna misalignment, and the correlation loss. The an- 
tenna loss is main ly related to our Cassegrain design. 
iKoch et al.l ()2009c( ) calculate the overall antenna effi- 
ciency to be 0.58 for the 1.2 m and 60 cm dishes alike. 
The loss originates mainly from three factors: the illu- 
mination efficiency, the secondary blockage, and the for- 
ward spillover. The AMiBA feeds provide a Gaussian 
illumination pattern and the the reflectors are designed 
to have a -10.5 dB edge taper. Compared to an uniformly 
illuminated reflector, only 90% of the dish is effectively 
used. The shadow of the secondary mirror blocks about 
8% of the collecting area, giving a loss of 0.92. The edge 
of the secondary corresponds to the edge of the primary 
mirror. Therefore the illumination is either reflected by 
the primary to the sky or is emitted toward the sky dire- 
cly. The latter part constitutes about 22% of the energy 
giving the forward spillover factor of 0.78. We favored 
a design with slightly worse forward spillover but lit- 
tle to none backward spillover (illumination toward the 
ground) to reduce system temperature. 

The antenna misalignment consists of the mechanical 
installation error and the dynamical deformation of the 
platform. The former error was measured to be around 
3' during the 2007 observing season (Wu et al. 2009, 
in preparation) and will be improved for future observa- 
tions. The latter error was inferred from photogramme- 
try measurements of t he platform surface to be less than 
1' (|Koch et al.ll2009al) . The two errors together atten- 
uate the primary beam by 2%. Antenna misalignment 
may also cause pointing errors for some baselines. This 
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Fig. 3. — The upper two rows display three complex response functions. From the left, two slightly different functions are taken from the 
measurements shown in Fig. [T]with all four lags plotted explicitly. The third is an assumed flat response function. We use the left-most 
function to simulate a set of fringes from a flat spectrum source moving along the baseline. Then we follow the lag-to-visibility procedures 
in i|2.4l to obtain the bottom three sets of complex visibilities using each of the response functions as the transformation kernel. The 
horizontal axis is the source offset presented as drift time. To the left we see the case when the kernel is the same as the response function. 
The central plots show the results with small errors in the kernel. On the right we show the result of using an assumed flat kernel. Even 
though the deviation from ideal behavior can be as large as 50% in visibilities, calibration can correct for the errors. 



effect is not considered in individual baseline efficiencies 
but will be considered in the array efficiency. There is 
also a loss of efficiency from the noise contributed by 
the rejected correlations in the analog correlator. The 
estimated correlation efficiency from this effect is 0.8f . 

When combining baselines from the entire array, point- 
ing error and system stability also lower the efficiency by 
degrading the coherence of signal from differe nt measure- 
ments . The pointing error is less than 0.4' ()Koch et al.l 
l2009al) and decreases the efficiency by less than 2%. The 
large alignment error, on the other hand, contributes as 
much as 12% loss in the 2007 and 2008 observations. As 
described in i j2.5[ the system stability is approximately 
±5% in gain and ±0.1 rad in phase. Taken over all base- 
lines this results in a reduction of signal by about 2%. 
Table [T] summarizes the major losses in the system. 

The baseline efficiency has been checked by comparing 
the signal-to-noise ratio (SNR) of Jupiter's fringe to the 
ratio of Jupiter's antenna temperature and the system 

temperature, -qu — SNRj„p/(2|t^ WSg/z^X^), where 



TABLE 1 
Summary of Losses of The System 



Systematics 


Efficiency 


Antenna illumination^ 


0.90 


Antenna blockage^ 


0.92 


Antenna spillover'' 


0.78 


Antenna others effects^''' 


0.90 


Antenna total 


0.90 X 0.92 X 0.78 x 0.90 = 0.58 


Alignment'' 


0.98 


Correlation'' 


0.81 


Overall Baseline 


0.58 X 0.98 X 0.81 = 0.46 


Deformation/Pointing'' 


0.88 


Stability'' 


0.98 


Overall Array 


0.46 X 0.88 X 0.98 = 0.40 



^ The facor is based on theoretical calculation.'' The facor is 



derived from measured quantities. 

Ta,jup is the antenna temperature of Jupiter, typically 
around O.IK for the 60cm dishes, and SNRj„p is the SNR 
of Jupiter under the corresponding recording time tree 
(= 0.452 sec currently). An average effective bandwidth 
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Fig. 4. — Visibilities recovered from a repeated two-patch tracking of Saturn (open triangles and plus symbols) and Jupiter (open circles 
and cross symbols) in the local time range 8 pm to 8 am (UT 6hr to 18hr). The upper panel shows the relative gain fluctuation, and 
the lower panel shows the phase variation. Both are shown for the XX correlations. The number in the upper-left corner of each subplot 
indicates the antenna combination. The scale of the plot and the UT time range are indicated at the bottom left subplot of each panel. 

of Be// ^ 10 GHz was assumed in the calculation. The 
measured efficiency scatters from 0.2 to 0.5 with an error 
bar of approximately 0.2. The error originates mainly 
from the noise estimation of the signal-dominant fringe, 
the occasional large readout noise, and also the variation 
of effective bandwidth. The overall array efficiency will 
be covered in 

3.2. Calibrator 

The raw visibilities recovered from the lag data have 
the systematic losses discussed above and are further af- 



fected by instrumental delay, gain drift, phase variation 
as well as the imperfect lag-to-visibility transformation. 
We calibrate visibilities by interspersed two-patch obser- 
vations of a planet. This converts our visibility ampli- 
tudes to flux density units and references the phase to 
the calibrating planet position. 

Taking planet data with the subtraction scheme, and 
applying the same calibration scheme used for cluster 
data (one calibration about every three hours), we find 
that the recovered peak flux in the image domain shows 
an RMS scatter of about 3%. 
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Fig. 5. — Jupiter flux recovered from the stability measurements 
in Fig|4] Two horizontal dotted lines indicate ±4% of the expected 
Jupiter flux. 

The flux densities of the planets are calculated from 
published disk brightness temperatures and the appar- 
ent angular sizes assuming a black-bo dy spectruin. We 
adopt the va l ues: J upiter 171.8±1.7K fiPage et al."2003'; 
iGriffin et al.][l986l ). Sat urn 149 3±4 1 K (Ulich 1981), 
and Mars 206.8±5.6 K (jUlich I Il98l . Fig [6] shows the 
recovered flux of the main calibrators (Jupiter and Sat- 
urn) in the 2007 and 2008 observations and the expected 
flux from calculation. For the phase reference and flux 
standard we use the first Jupiter measurement on each 
night when Jupiter is observed, as this provides the best 
calibration for the earlier Saturn measurements. The 
scatter of the Jupiter flux densities agrees with the scat- 
ter in one night as shown in FigO The flux density of 
Saturn is systematically lower than the calculated value 
by approximately 5%. It was verified that the lower flux 
was not the result of an error in flux standard. This 
was checked by artificially setting the phase error to zero 
in the calibrated visibilities of Saturn and forming an 
image. The recovered flux density displays no system- 
atic offset from the calculated level, to within the error 
bound. The apparent deficit in the fiux density for Sat- 
urn is likely due to insufficient phase calibration in the 
2007 data, when calibrations more than two hours apart 
were often used. Note that the effect of Saturn's ring is 
not included in the calculation of Saturn's flux density. 
The ring inclination was about -15 dog for the 2007 ob- 
servations and about -10 deg for the 2008 observations. 
We estimate that our flux density scale is good to about 
±5% in absolute terms. Based on Jupiter's flux scale, our 
measurement of Saturn's disk brightness temperature is 
149+^2 K. 

A final note about calibration is on the difference of 
spectra between the calibrator and the SZ effect. Our 
primary calibrators, Jupiter and Saturn, are dominated 
by thermal emission near 94 GHz. The black-body spec- 
trum favors slightly higher frequency than does the SZE 
spectrum. Regardless of the choice of assumed passband 
shape, the difference of effective central frequency, de- 
fined by /, = i?(/)5(/)/d/]/[/,^„, i?(/)5(/)d/] 

with S{f) being the source spectrum and R{f) being 
the passband, is less than 1.5% of the bandwidth. Its 
effect on the calibration is thus negligible compared to 
other errors. 

3.3. Noise Integration 
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Fig. 6. — The recovered flux densities of our main calibrators, 
Jupiter (plus symbols) and Saturn (cross symbols) from the obser- 
vations in 2007 and 2008. Data were calibrated by the first Jupiter 
measurement in the same night, or the nearest Jupiter measure- 
ment. The expected flux densities are plotted as dashed and solid 
linos for Jupiter and Saturn respectively. Note that Saturn's ring 
was not taken into account in the flux density estimation. 

Based on the signal-to-noise ratio (SNR) of Jupiter's 
fringe and the discussion in ^JSHJ we find that the ar- 
ray has an overall efficiency of about 0.4. The im- 
provement in sensitivity with integration time depends 
critically on the removal of spurious signals using the 
subtraction scheme and subsequent data flagging. To 
verify the sensitivity, we examine the variation of sig- 
nal and noise in the reconstructed map with integration 
time in Figure [7l The integration time here refers to 
the accumulation of observing time spent in each in- 
dividual visibility band. For example, when the tele- 
scope tracks a source for 3 minutes, the total integra- 
tion time for 21 baselines, 2 polarizations, and 2 bands 
is ttot = 180 sec X 21 X 2 X 2 = 15120 sec. Since the 
visibilities are used with non-uniform weights in form- 
ing the image, we calculate an effective integration time, 
which for 3 minutes on-source integration is deflned as 

tf,ff = ^1^' X 180 sec, where Wi denotes the weight- 



XT^ 



ing given to each data set. In this analysis, a natural 
weighting is adopted. 

Since the noise comes from two patches and the signal 
comes from only one, the point source sensitivity can 
be estimated by cr = , ^ where T,„, = 

100 K, Bch = 5 GHz, T]aii is the overall efficiency, and 
Aphys is the physical collecting area. 

The signals are read from the source position in the re- 
constructed dirty images with different integration times, 
with the source p osition determin ed from the flnal im- 
age. A CLEAN (jHogbom Ill974[) procedure is applied 
to the inner 21.6' box, which roughly corresponds to the 
FWHM of the primary beam. The cleaned signal at the 
source position is also recorded, while the residual noise 
is measured in the 1 deg image excluding the inner clean 
region. 

Figure [7] shows that an efficiency rj = (0.36 ± 0.04) 
is more representative for the current data, giving a 
point source sensitivity of 63 ± 7 mJy in 1 hour of on- 
source integration when the subtraction scheme is ap- 
plied. Abell 2390 has a higher noise, and we believe this 
may be caused by th e presence of point sources in the 
fov. iLiu et all (|2009( ) investigate the contamination by 
point sources and the primary CMB. 



8 



Lin et al. 



E 
E 



1000 



100 



1000 



w 

CD 

B 100 

.D3 



1000 



100 



A1689 



□ □ 




A2142 




A2261 



A1995 



peak(dirty) □ 

peak (cleaned) ■ — ^ 

noise (1 a) 

noise (ri=0.36) 




A2163 




A2390 




1000 10000 100000 1e+06 1000 10000 100000 1e+06 

integration (sec) 
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4. CONCLUSION 

To detect galaxy clusters with the AMiBA, we must 
achieve system stability on timescales of hours. We have 
optimized the performance of AMiBA by measuring and 
compensating for the delays between antennas, and using 
beam switching techniques to cancel out instrumental 
and environmental effects. Planet calibrations provided 
corrections for passband response. Overall efficiency for 
AMiBA was ijaii — 0.36 ± 0.04, with a major loss from 
the antenna efficiency, rjant = 0.58. 

Using a system temperature of 90 K, an effective band- 
width of 5 GHz per band, and an overall efficiency of 
0.36 ± 0.04, the point source sensitivity of AMiBA in 
1 hour of on-source integration {teff = 302400 sec) is 
found to be about 63 ± 7 mJy when the subtraction 
scheme is applied. The effective integration is about 60% 
of the on-source integration time. The loss of 40% of ob- 
serving time is mostly due to lower weighting applied to 
some receivers or baselines experiencing hardware prob- 
lems. There were very few observations made when the 
weather was not good in the 2007 and 2008 observing 



season. 

The flux density error consists of the calibrator flux 
scale uncertainty of ±5 %, and the cross-calibration error 
±3%, which also includes the lag-to-visibility flux scale 
error of ±2%. The latter two errors are well below the 
thermal noise in all clusters observed during 2007 and 
2008. Investigation of noise in cleaned images shows that 
longer integration, aimed at measuring primordial CMB 
fluctuation, will be promising. 

The resulting successful detections of clusters have led 
to a number of scientiffc results including a measurement 
of the Hubble constant and the study of the hot gas dis- 
tribution in the clusters. These are discussed further in 
the companion papers. 



Acknowledgments Capital and operational funding 
for AMiBA came from the MoE and the NSC as part of 
the Cosmology and Particle Astrophysics (CosPA) ini- 
tiative. Additional funding also came in the form of an 
Academia Sinica Key Project. Support from the STFC 



AMiBA Performance 



9 



for MB is also acknowledged. 



Bond, J. R. et al. 2005, ApJ, 626, 12 
Chen, M.-T. et al. 2009, ApJ, submitted 
Goldstein, J. H. et al. 2003, ApJ, 599, 773 

Griffin. M. J. ct al. 1986, Icarus, 65, 244 

Ho, P. T. P ct al. 2009, ApJ, accepted {astro-ph/0810.1871) 

Hogbom, J. A. 1974, A&AS, 15, 417 

Huang, C.-W. L. 2009, ApJ, submitted 

Huang, Y.-D. ct al. 2008, SPIE, 7012, 70122H 

Koch, P. M. ct al. 2009a, ApJ accepted 

Koch, P. M. ct al. 2009b, ApJ, submitted 

Koch, P. M. et al. 2009c, in preparation 

Kuo, C.-L. et al. 2004, ApJ, 600, 32 

Li, C.-T. et al. 2004, SPIE, 5498, 455 



REFERENGES 

Lin, K.-Y., Woo, Tak-Pong, Tseng, Yao-Huan, Lin, Lihwai, and 

Ghiueh, Tzihong 2004, ApJ, 608, IL 
Lin, K.-Y. et al. 2008, SPIE, 7012, 701207 
Liu, G.-C. et al. 2009, ApJ, submitted 

Nishioka, H. et al. 2009, ApJ accepted {astro-ph/0811.1675) 

Padin, S. et al. 2002, PASP, 114, 83 

Page, L. et al. 2003, ApJS, 148, 39 

Pearson, T. J. ct al. 2003, ApJ, 591, 556 

Umctsu, K. ct al. 2009, ApJ accepted (astro-ph/0810.0969) 

Wu, J.-H. P. et al. 2009, ApJ accepted (astro-ph/0810.1015) 

Ulich, B. L. 1981, AJ, 86, 1619 



